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We study plasma effects in a micromachined high-electron mobility transistor (HEMT) with the 
microcantilever serving as the gate using the developed a model. The model accounts for mechanical 
motion of the microcantilever and spatio-temporal variations (plasma effects) of the two-dimensional 
electron gas(2DEG) system in the transistor channel. The microcantilever mechanical motion is 
described in the point-mass approximation. The hydrodynamic electron transport model is used to 
describe distributed electron plasma phenomena in the 2DEG system. Using the developed model, 
we calculated the response function characterizing the amplitude microcantilever oscillations and 
the output electric signal as functions of the signal frequency and the bias voltage for the devices 
with different parameters. We find the voltage dependences of the frequency of the mechanical 
resonance and its damping. In particular, it is demonstrated that the amplitudes of the mechanical 
oscillations and output electric signal exhibit pronounced maxima at the bias voltages close to the 
voltage of the 2DEG channel depletion followed by a steep drop with further increase in the bias 
voltage. 

PACS numbers: 73.50. Mx, 73.40.-c, 73.43.Cd 
PACS numbers: 



I. INTRODUCTION 

The concept of a field-effect transistor with a micro- 
machined cantilever as floating gate was put forward 
and discussed a long time ago by Nathanson et al [l| 
(see, also ref. 0). Recently 0, this concept was fa- 
ther evolved and a floating-gate high-electron-mobility 
transistor (HEMT) device comprising a microcantilever 
over a two-dimensional electron gas (2DEG) channel 
was fabricated and measured (see, for instance, also 
refs. [1, 0, [E IE A theoretical assessment of this 
and similar devices is usually based on equations gov- 
erning the mechanical motion of the cantilever supple- 
mented by some circuit equations. However, the mo- 
tion of the electrically actuated microcantilever results 
in a complex motion of electrons in the 2DEG channel, 
which, in turn, leads to variations of the electron sheet 
density and the self-consistent electric field affecting the 
microcantilever dynamics. The simplest effect of 2DEG 
in micromachined HEMTs is associated with a finite con- 
ductivity of the latter, particularly at the voltages corre- 
sponding to a significant depletion of the 2DEG channel. 
The electron transient processes associated with the vari- 
ations of the electron sheet density and the self-consistent 
electric field, i.e., electron plasma processes, can play a 



significant role in HEMTs and different HEMT-based de- 
Il3ll (see also experimental pa- 
QjJ 120] )■ The plasma effects in 



12 
18 



'Electronic address: v-ryzhii@u-aizu.ac.jp 



the HEMT-like devices can result in resonant response at 
the frequencies coinciding with the frequencies of plasma 
oscillations. Since these frequencies usually fall into the 
terahertz range, such effects are important in different 
terahertz devices. Despite a substantial difference in the 
resonant mechanical frequencies and the characteristic 
plasma frequencies, micromachined elements such as a 
microcantilever floating-gate may be useful for terahertz 
devices. For instance, the microcantilever floating gate 
can be used for the internal modulation of terahertz sig- 
nals in the resonant detectors utilizing the excitation of 
plasma oscillations (see, for example, ref. [8] ) . The con- 
sideration of plasma effects in the HEMT-like devices re- 
quires a device model which adequately describes spatio- 
temporal variations of the 2DEG system. 

In this work, we develop a model which self- 
consistently describes the mechanical oscillations of a 
highly conducting (metallized) microcantilever and the 
dynamic properties of the 2DEG system. The model 
is based on an equation governing oscillations of a mi- 
crocantilever under mechanical and electric forces and 
hydrodynamic equations governing the electron trans- 
port in 2DEG. Our model accounts for such phenom- 
ena as the depletion and enrichment of 2DEG by the 
applied voltage in the presence of the surface charges at 
the semiconductor surface beneath the microcantilever 
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and finiteness of the 2DEG conductivity, the delay in 
the electron recharging under the microcantilever (gate) 
and the spatial nonuniformity of the potential distribu- 
tion in the 2DEG, which might affect the device char- 
acteristics. Thus, the model under consideration is a 
distributed (physical) model which provides more gen- 
eral and detailed description of the underlying processes 
than that based on the treatment of the electron system 
of the device invoking electric circuit models. Generally, 
the device under consideration is similar to those fabri- 
cated and studied both theoretically and experimentally, 
in particular, in r efs. [HQ. The main distinction is that 
in micromachined HEMTs investigated in refs. [l|,l2] there 
are two metallized sufraces: one at the bottom of the mi- 
crocantilever and one covering a portion of the semicon- 
ductor surface under the microcantilever (the so-called 
input force plate [l[). In contrast, we assume that the 
microcantilever motion is due to the interaction of the 
charge induced in the microcantilever metallized surface 
and the charge in the 2DEG channel. Thus, the 2DEG 
channel plays the dual role: it is used for the micro- 
cantilever actuation and the variation of its conductivity 
is used to detect the output signals. The effectiveness 
of the control of the 2DEG channel by applying voltage 
to the microcantilever was demonstrated experimentally 
even in ref. [l[ (see, also ref. 0] and references therein). 
The effect of interaction of the charges induced in the 
microcantilever and the conducting plane electrode was 
theoretically considered recently in ref. 21]. However, 
the model used in ref. [21] assumes, in contrast to our 
model, that the electrode is ideally conducting that is 
not always the case in real 2DEG channels. 

The developed model is used to study the resonant 
response of the device to the ac signals. We find the de- 
pendences of the resonance frequency and the resonance 
width on the both mechanical and electron properties of 
the system under consideration. In particular, we demon- 
strate that the plasma effects in the 2DEG channel can 
give rise to a shift of the resonance and an increase in the 
oscillation damping, i.e., limit the quality factor of the 
microcantilever oscillations. We show also that relatively 
high-frequency signals (with the frequency corresponding 
to the plasma resonance) can result in a significant vari- 
ation of the microcantilever position. 

The paper is organized as follows. In Sec. 2, we discuss 
the device model and write down the pertinent equations. 
Section 3 deals with a small-signal analysis based on the 
linearized versions of equations of the model. In this sec- 
tion, we calculate the spatial distributions of the ac po- 
tential and electron density in the 2DEG channel and find 
the amplitude (response function) of the microcantilever 
oscillations as a function of the signal frequency and the 
bias voltage. In Sec. 4, we analyze the response function 
in different limiting cases and demonstrate the results of 
numerical calculations of the device characteristics using 
the formulas obtained in the previous sections. Section 5 
deals with the calculations of the output source-to-drain 
ac current and the effective device transconductance. In 
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FIG. 1: Schematic view of device structure. 
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FIG. 2: Device model. 



Sec. 6, we consider the effect of high-frequency signals ( 
on the variation of the microcantilever position associ- 
ated with the plasma resonance. In conclusion (Sec. 7), 
we draw the main results. Some ancillary calculations 
are factored out to Appendix A and Appendix B. 



II. EQUATIONS OF THE MODEL 

The device structure under consideration is schemati- 
cally shown in Fig. 1. It is assumed that the voltage ap- 
plied between the metallized cantilever, which serves as 
the HEMT floating gate, and the side contacts to 2DEG 
channel (HEMT's source and drain) comprises the dc 
(Vo) and ac (SV) components. Focusing on fairly de- 
tailed description of the electron transport in the 2DEG 
channel accompanying the microcantilever oscillations, 
we shall consider the cantilever mechanical properties in 
the framework of a simplified model, namely, using the 
so-called point-mass model [12] (see Fig. 2). This model 
assumes that consideration of the elastic microcantilever 
or beam is replaced by the consideration of a point mass 
M (cantilever effective mass) attached to a string with 
stiffness K with M and K chosen such that the resonant 
frequency of the microcantilever oscillations associated 
solely with its mechanical properties, fio = \J K/M. In 
this model, the displacement of the cantilever (gate) is 



governed by the following equation: 
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dx£(E— Sd+S s ), 
(1) 



where Z = Z(t) is the distance between the cantilever 
surface and the 2DEG, W is this distance in the absence 
of the applied voltage in equilibrium, 70 is the damp- 
ing of the cantilever oscillations associated with different 
mechanisms of the energy loss in the cantilever body and 
in the clamp, e = \e\ is the value of the electron charge, 
D and L g are the pertinent sizes of the cantilever (see, 
Fig.2), £ = S(t, x) is the electron sheet density of 2DEG, 
£d = const is the donor sheet density, and E s = q s /e, 
where q s is the sheet density of the surface charge at the 
interface between the semiconductor and the gas (or vac- 
uum), which does not change when the microcantilever 
moves. The surface charge corresponds to the surface po- 
tential V s = 4:7rewT, s /se. The electric field £ = £(t,x) at 
the microcantilever plane is determined by the potential 
drop Vq — V s — <p, where ip — ip(t,x) is the potential of 
2DEG. Here, the axis z is directed perpendicular to the 
2DEG plane, while the axis x is directed in the 2DEG 
plane. Equation (1) is valid in the case of a highly con- 
ducting (metallized) cantilever when its surface is equipo- 
tential. The microcantilever and 2DEG are separated by 
two layers: the isolation solid layer of thickness w and 
dielectric constant as and the layer of a gas (or vacuum) 
of thickness Z — w and dielectric constant as* ~ 1 (see 
Figs. 1 and 2). The term in the right-hand side of eq. (1) 
represents the electric force acting on the microcantilever 
due the applied voltage. 

In the gradual channel approximation valid if Z, W <C 
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[Z — w(l — ae*/ae)] ae* 



£-£ d + £ s ). (2) 



Equations (1) and (2) disregard the effect of fringing ca- 
pacitance, [231 ] which can appear at the bias voltages be- 
yond the essential depletion of the gated portion of the 
2DEG channel (this voltage range is not considered here) . 
Equations (1) and (2) lead to 
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dx£ 2 , (3) 



where Z > w and £ > 0. In the following, we put for 
simplicity 88* = 1. 

Equations (2) and (3) should be supplemented by the 
hydrodynamic equations (continuity equation and Eulcr 
equation) governing the electron transport in 2DEG 0, 
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dt dx 



(4) 



Here, u = u(t, x) is the average (hydrodynamic) electron 
velocity in the 2DEG plane, v is the electron collision fre- 
quency, and m is the electron effective mass. The electron 
collision frequency is expressed via the 2DEG mobility fi 
as v — e/mfi. 



III. SMALL-SIGNAL ANALYSIS 

Let us assume that the net voltage between the can- 
tilever (gate) and the side contacts to the 2DEG (source 
and drain) apart from dc components (Vq and Vo + Vd) 
comprises also the ac component SV exp(— iuot): where 
SV (\SV\ <C Vq) and u> are the amplitude and frequency 
of the ac voltage. The drain-to-source voltage Vn is as- 
sumed to be sufficiently small and corresponds to the 
linear region of the HEMT operation. Considering small 
oscillations of the cantilever and the electron density, one 
can assume 

Z(t) — Zq + SZ exp(— iwt), 
£(i, x) = So + <5£ exp(— iu>t), 
u(t,x) = 5u exp(-itut), 

ip(t,x) = Sip exp(— iuit), 

where the amplitudes SZ, ST,, 5u, and Sip are assumed 
to be small in comparison with the steady-state sepa- 
ration, Zq, between the microcantilever and the 2DEG 
channel (see Appendix A). Here So = £d — T s + [(Vo — 
V s )/AireZ ] ~ E d - S s + (y /47reZ ) is the dc elec- 
tron density in the gated portion of the 2DEG chan- 
nel. In the above equation (and in the following), we 
have omitted for brevity the term with V s because of its 
smallness (it is proportional to small values w and ae -1 : 
V s oc wT, s /se). Nevertheless, the direct contribution of 
the surface charges to the dc electron density is taken 
into account. Considering the smallness of the above- 
mentioned amplitudes of variations, eqs. (1) - (5) can be 
linearized. As a result, neglecting w(l — ae -1 ) in compar- 
ison with Zq, we arrive at 
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dx 2 47re 2 So^o 47re 2 So^o Zq 

or introducing the characteristic plasma velocity 5* as 
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we obtain the following equations 



fi — fi ( = ) — ijoUJ — uj 



SZ 



L ^ 2 , 5<p 
dx——, 



Lg/2 



V 



and 



d 2 Sip uj(uj + iv) uj(uj + iv)Vo SZ 

-Sip = — 



S 2 



S 2 



(13) 



The boundary conditions for eq. (13) can be chosen 
considering that the amplitude of the ac potential of the 
side contacts is equal to SV and taking into account the 
ac potential drop across the access sections of the channel 
(see Fig. 2), i.e., the regions between the gate edges and 
the side contacts. Assuming that the electron collision 
frequencies in different regions of the 2DEG channel are 
the same and neglecting the deviation of the electron 
sheet density in the access regions from the donor density 
Srf, the boundary conditions can be presented as [24| 
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x=±L g /2 



Consider first the device structure in which the length, 
L a , of the access (ungated) sections, is sufficiently small 
in comparison with the gate length L g (i.e., L a <C 
LgT,d/T, ). In this case, one can neglect the ac potential 
drop across the access regions of the 2DEG channel, i.e., 
the second term in the right-hand side of eq. (14) fillip 
As a result, from eqs. (13) and (14) we obtain 
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Substituting Sp from eq. (15) to eq. (12), we arrive at 
5V = V Q Z "' (16) 

where 

Z^ = 



nKVo/VofiUnQ^/Q^) 



[fig - i l0 uj -uj 2 - fig(^o/^o) 2 (tang w /Q w )] 



Here, we have introduced 
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where 



TT \J Uj(bJ + iv) 



v a = 




^fi 2 MW^ (18) 



(19) 



(12) is the characteristic plasma frequency of the gated 2DEG 
channel and Q p0 = y / 47r 3 e 2 (E d — T, s )W/mL 2 

One can see that fi p depends on the bias voltage via 
the voltage dependence of £o and Zq. The fi p versus Vq 
dependence is asymmetric; The plasma frequency can be 
significantly decreased by negative bias when the 2DEG 
channel becomes close to the depletion. 

The finiteness of the conductivity of the access ungated 
regions can also contribute to the damping of the micro- 
cantilever oscillations. This can occur when the length 
of this regions L a is sufficiently large. To include these 
regions into the model, we need to modify boundary con- 
dition (19) to take into account the potential drop across 
them. Generally, the access regions can pronouncedly af- 
fect the plasma phenomena in HEMTs [HI, 0, , in par- 
ticular, leading to a decrease in the characteristic plasma 
frequencies. 

Preserving the second term in the right-hand side of 
eq. (4) associated with the contribution of the access re- 
gion to the boundary conditions, we arrive at 



n 2 (v /Vo) 2 {^Q^/Q*J 



{fi 2 - i lQ uj -uj 2 - fi 2 (Vyy ) 2 (tanQ w /Q*)} 

(20) 

where Q* = Q w (l — aQ^tanQ^) and a = 
(2L a / L g )[Y,o/(T, c [ — E s )] is the parameter characterizing 
the role of the access regions. At a = 0, Q* = Q u and 
eqs. (17) and (20) coincide. 



IV. MICROCANTILEVER FORCED 
OSCILLATIONS (ANALYSIS OF LIMITING 
CASES AND NUMERICAL CALCULATIONS) 

A. Highly conducting 2DEG channel 

In many practical situations, the signal frequency uj is 
in the the same range as the resonant frequency of the 
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microcantilever oscillations f2o and the conductivity of 
the 2DEG channel is rather large. The latter corresponds 
to Q p > v. Since usually u>, fio ^ ^jQp, the quantity 
\Q U \ <C 1, and eqs. (17) and (20) can be simplified. In 



such a case, one obtains t&nQ u /Q u 



1 



1 



i(ir 2 /l2)(ujv)/flp). Considering this, assuming that a <C 
1 and using eq. (17), the response function Z w can be 
presented in the standard form 



n 2 

s 'o 



(21) 



Here, however, the resonant frequency O and the quantity 
characterizing the damping of oscillations 7 depend on 
the "electron" parameters: 



7m = 70 




(22) 
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(23) 



where (3 = Vb/|V r Q dep/ ' l | and y o (acpl > i s the characteristic 
depletion voltage or HEMT's threshold voltage. Here 
we have considered the dependence of the characteristic 
plasma frequency on S and the dependence of the latter 
on the bias voltage (see Appendix B, cq. (B2)). As seen 
from eq. (21), the resonant frequency is equal to Q. The 
modulus of the response function \ZJ\ is given by 



(depl) 



o 2 
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so that at the resonance 
max IZJ = 



uj 2 ) 2 +jIuj 2 Wr 



Q 2 
s 'o 



7m fir. 



(24) 



(25) 



The second term in the right-hand side of eq. (23) pro- 
portional to v/Q 2 oc f/So (i- e -> proportional to the re- 
sistance of the gated region of the 2DEG channel) deter- 
mines the contribution to the resonance width associated 
with the dissipation processes in the gated 2DEG channel 
due to the finiteness of its conductivity. Equation (24) 
demonstrates a tendency for 7m to increase with ap- 
proaching to the 2DEG channel depletion. If the second 
term in the right-hand side of eq. (23) becomes dominant, 
it might limit the quality factor of the microcantilever os- 
cillations: Q m = fl m /j m < Qo = £^0/70 • F° r example, 
for n /2n = 100 MHz, Q p / 27r = 50 GHz, v = 4x 10 12 s' 1 
(electron mobility fi = 8000 cm 2 /Vs), and V /% = 0.1, 
that can correspond to a HEMT with GaAs channel close 
to depletion at room temperature, the quality factor is 
limited by the value maxQ m < 5 x 10 3 . 



B. Low conductivity of the 2DEG channel (low 
plasma frequency) 

At a strong depletion of 2DEG channel at negative 
bias voltages, the conductivity of 2DEG channel and the 
plasma frequency can be relatively low. At low plasma 
frequency fl p when lu,Qo 3> Q 2 ,/^, one obtains ~ 
7r\/ iuuv/2VL p with \Q U \ = TTy/uJi'/2ilp 3> 1. In such a 
situation, tanQ^/Q^ ~ [(1 + i)\/2/n](i} p /^/Uuj). Hence, 

z n 2 (V /Vo) 2 [(l + i)V2n p /TrVl^] 

" ~ n 2 - i loJ - lu 2 - n 2 (v /"^) 2 [(i + i)V2~n p /ir^J] ' 

(26) 

As follows from eq. (26), when Q ^ ^>l v at the reso- 
nance (compare with the pertinent formula for the case 
highly conducting 2DEG channel), 



maxlZ^ 



fio 

7o 



2 n 



v 



fio 

7o 



(27) 



Thus, in the case under consideration here, even at 
the resonance, max|Z w | < 1. If y/ u/Hq, increases, 
i.e., the 2DEG channel conductivity decreases, max|Z w | 
markedly decreases as well. Hence, when the bias volt- 



age Vo approaches to the depletion voltage Vq 
resonant peak markedly diminishes. 
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C. Mechanical response at plasma resonance 

If the signal and plasma frequencies markedly ex- 
ceed the electron collision frequency and the resonant 
frequency of pure mechanical oscillations of the micro- 
cantilever (w, Q p 3> v, f2o) and the signal frequency is 
close to one of the plasma frequencies Q p (2n — 1), where 
n = 1,2, 3,... is the index of the plasma mode, the quan- 
tity |tanQ w /Q w | can be rather large. At the fundamen- 
tal plasma resonance tanQ^/Q^ ~ i(4/ir 2 )(fl p /is). Tak- 
ing this into account, from Eq. (18) we find that as in the 
case considered in the previous subsection, max|Z w | < 1. 
One can also find that when Q, p /v increases to infinity, 
maxlZ^I formally tends to unity. However, in most real- 
istic cases (Qo ^<fi P and Vq < Vq), 



max|Z w | ~ 



n 2 



< 1. 



(28) 



This shows that even at rather sharp plasma resonance, 
the amplitude of the microcantilever oscillations remains 
small in comparison with the amplitude at the mechani- 
cal resonance at u = Cl. Indeed, the ratio of max|Z w | at 



lu = Q p and at u) = O, can be estimated as 



maxIZ. 



=fi p 7T 2 7 Q 



maxIZJ^n 4 



< 1. 



(29) 



At the plasma resonance, the amplitude of the ac electric 
field in the gated region can be rather large. 
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D. Role of the access regions 

Consider now the case when a is not small invoking 
Eq. (20). When u),Qq <C u, Cl p , one obtains \QJ\ — 
■k^/ujv /2Q p <C 1 and eq. (20) can be simplified. In the 
vicinity of the resonance to ~ f2 m , where in the case uder 
consideration (compare with eq. (22)) 



TO — Qq * / 1 — 



(Vo/V 



[l + a 2 ^) 4 ^) 2 /^ 4 ]' 



one obtains 

Z, 



n 2 o(V Q /V ) 2 



(fl 2 - ijLU - u 2 )[l - ia(Tr 2 /4)(LOv)/n 2 p } 



(30) 



(31) 




Here 



FIG. 3: Modulus of the response function \Z^\ as a function 
of normalized signal frequency ui/SIq at different bias voltages 
0Vo/Vo = V o /\v£ depl) \. 



r 7r 2 (l + 3a) 



7m = 70 + ^ 
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(n /n p ) 2 fv Q 

[1 + ^(^/2)4(0^)2/04]^ 



7o + v 



7r 2 (l + 3a) 
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(n /n p ) 2 

(1 + ^r 2 ) 



(32) 



is the quantity characterizing the damping of oscillations. 
Here r a = a(ir 2 /4)v/S} 2 is the RC-delay time: r = 
R a C g , where R a — {mvL a /2e 2 Y>d) and C g — Lg/inZo 
are the resistance of the access regions and the effective 
capacitance of the gated portion of the 2DEG channel, re- 
spectively. When a(n 2 /4)(n v)/ill = ft T Q ~ 1, eq. (32) 
results in 



7o + v 



TT 2 a 



Oo 
{ 2 1, 



(33) 



One can see that the second term in the right-hand side 
of eq. (33) differs from the pertinent term in eq. (23) by 
a factor 3a/2 which can be large (when L a ^>> L g ). 
Equation (31) leads to 



O 2 
s 'o 



(34) 



At the exact mechanical resonance w = fi m , from 
eq. (34), we obtain 



maxlZ,, 



n 



V 



lm y/i + n 2 r 2 \v 



(35) 



One needs to point out that the second terms in the 
right-hand side of eqs. (24) and (34) are the products 
of small factors, (f^o/^p) 2 an d (Vq/Vq)2, and relatively 
large value v (normally v ^> 70). 



E. Results of numerical calculations 

Figures 3 shows the modulus of the response func- 
tion \ZJ\ determined by Eq. (20) versus signal frequency 
calculated for AlGaAs/GaAs micromachined HEMT de- 
vices under consideration at different bias voltages. The 
device parameters are chosen as follows: f2o/27r = 
100 MHz, Q = ^0/70 = 2500, fi p0 / 27r = 1 THz, and 
v = 10 12 s _1 . The above parameters correspond a Al- 



GaAs/GaAs HEMT-based device with to L g 



/mi, 



1 fim, W = 0.5 /mi, M/L g D = 4x 10" 12 g//mr 



1 x 10 cm , and the electron mobility /t 



s Q 

3x 10 4 cm 2 /Vs. At the above parameters, f2 p o/^o = 10 4 , 
u/Q,Q = 10 4 /2tt, % = 33 V, V$ depl) = -9 V (so that 

P = 3.67), and v£ pull ~ in) = 18 V. It is seen from Fig. 3 
that the resonant frequency decreases with increasing 
bias voltage (in line with eq. (22)). Figure 3 also shows 
that \Zu\ increases with increasing Vq reaching a maxi- 
mum at the voltage slightly smaller than the depletion 
voltage V ep . Further increase in Vq leads to a drastic 
drop in This is confirmed by Fig. 4 which shows the 
voltage dependence of max \ZJ\ calculated for the above 
parameters (curve 1). The results of calculations for a 



device with Sw 



10 1 



and L„ 



2 jj,m) 



are shown in Fig. 4 as well. The parameters related to 
curve 1 in Fig. 4 correspond to |v^ dep ^| < Vq PuL '"' ; " 



contrast to curve 2 for which \V r 



(depl) I 



(pull — in) 







V. OUTPUT ELECTRIC SIGNAL 

The ac voltage applied between the microcantilever 
and the source and drain contacts and the microcan- 
tilever oscillations both result in the modulation of the 
electron density in the 2DEG channel and, therefore, in 
the modulation of the source-to-drain current. Indeed, 
using eqs. (7) and (15), one can obtain (for the case of 
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FIG. 4: Maximum (resonant) value of \Z^\ as a function of 
normalized bias voltage /3 Vo/Vo = Vo/\Vq \ calculated for 
different structural parameters: 1 - E<j — E s = 1 x 10 11 cm -2 
and 2 - E d - E a = 3 x 10 11 cm" 2 . 



short access regions and relatively small signal frequen- 
cies) 



~ (X + Z u) cos[V^(tJ + iv) x/S] §v 
4vreZ cos[y/u}(uj + iv) L g /2S] 



The effective transconductance g m includes both the 
usual component associated with the direct electron den- 
sity modulation by the ac voltage and the component 
associated with the microcantilever oscillations (propor- 
tional to ZJ). In the limit of highly conducting 2DEG 
channel, eq. (39) can be presented as 



V d 



9n 



go |U(*=p0| 



1 



o 2 

"0 



V 



(fi 2 - i^/uj - uj 2 ) \ Vi 



(40) 



As seen from eq. (40), at the resonant frequency, the 
modulus of the transconductance exhibits a rather high 



maximum: 



max- 



Vp gg fVo 



go \V( de P l ) \jfl\VoJ \V^ 1 ) | 70 \V 



v D n rv c 



Since in the devices with high quality factor of mechan- 
ical resonance \Z U \ at the resonant frequency can be 
much larger than unity, the output signal in the microma- 
chined HEMT can significantly exceed the output signal 
in HEMTs with solely electrical modulation However, as 
the bias voltage is approached to the depletion voltage, 
the contribution of the microcantilever oscillations to the 
transconductance vanishes and the latter steeply drops 
to zero. As follows from eq. (41) and Figs. 3 and 4, the 
mechanical resonance provides very sharp peaks of \Z U \ 
and, hence, the transconductance modulus. 



(1 + ZJ) 



AneZo cos[ViuJi> L g /2S' 



5V ~- {1 + Zuj) 5V. 



AneZo 



At the signal frequencies close to the microcantilever res- 
onant frequency, i.e., at the frequencies significantly lower 
than the HEMT characteristic frequencies, one can use 
the transistor steady-state characteristics. As a result, 
the ac component of the drain current, which can be 
considered as the output signal, at low drain-to-source 
voltages Vd -C \Vq\ is given by 



6J d = e -^6X. 



Using eqs. (36) and (37), one obtains 

9oV d M , ~ x „, 9oV d 



5J d = 



AneYsdZo 



(1 + Z W )SV 



y(depl) 



(37) 



(38) 

where go — efiSd/Lg is the conductance of the unde- 
pleted 2DEG channel. Accounting for eq. (17), from 
eq. (38) we obtain the following formula for the effec- 
tive transconductance g m = — (dJd/dV)\v d =const of the 
micromachined HEMT under consideration: 



9m 

9o 



V d 



y(depl) 



(1 + ZJ) 



V d 



|y(dep0| 



(Vo/Vo) 2 (tan Q UJ /Q a 



VI. ELECTRIC-FIELD OSCILLATIONS NEAR 
THE PLASMA RESONANCE 

Using eqs. (2) and (15), one can find the ac electric 
field 88 as a function of 8V and the signal frequency oj. 
As shown the previous section, in the range of signal fre- 
quencies where the ac electric field can exhibit the plasma 
resonances (u> ~ Sl p S> Oo, v), the ac displacement of the 
microcantilever is relatively small. Hence, one can ne- 
glect the first term in the right-hand side of eq. (16) and, 
as follows from eq. (2), put 8£ ~ Scp/Zi. Here, Z\ = (Z) 
is the dc position of the microcantilever under the dc elec- 
tric field between the microcantilever and 2DEG channel 
and the average effect of the ac electric field (compare 
with eq. (A2)): 



Z 1 = W 



DL„ 



Antt 2 M 



V 



1 



dx(£ 2 



-L g /2 



(42) 



where the symbol (...) means averaging over fast oscilla- 
tions with the frequency to 3> ^o- As a result, we arrive 
at the following equation: 



5£ 



cos[-y/w(w + iv) x/S] SV 
cos[-\/cj(u; + iv) L g /2S] Z\ 



(43) 



[fig - i-iooj -oj 2 - 9,l{Vo/Vo) 2 {t&nQu/Qu>)} 

(39) 



Substituting S£ from eq. (43) into eq. (42), averaging 
■over high-frequency oscillations, we obtain 
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Here 



4^ J-L a /2 



cos[^/w(w + iv) x/S] 



cos[-\/w(w + iv) L g /2S] 



(45) 



At relatively low dc and ac voltages when (W — 
Z )/W < 1 and (W - Z x )/W < 1 (see cq. (A4)), we 
obtain the following formula: 



W 



VoJ \VJ 



(46) 



The second and third terms in the right-hand side of 
eq. (46) are associated with the microcantilever displace- 
ment due to the bias voltage and the ac signal, respec- 
tively. Assuming Q p 3> v, we obtain 



[1 + {2n p /iruj) sin(7ru>/20 p ) cos(ttu;/20 p )] 
8[cos 2 (7Tw/2rj p ) + (irv/4il p ) 2 } 



(47) 



As follows from eq. (47), the electric response function 
T exhibits sharp resonances at u — fl p (2n — 1), where 
n = 1,2,3, ... is the resonance index, if Q p 3> v. In the 
vicinity of the fundamental resonant frequency, i.e., at 
u ~ O p , eq. (41) yields 



1 



9 2 



2tt 2 [(oj-n p ) 2 + v 2 /A\ 
At this resonance, one obtains 

2 /Q, 

maxJV, ~ — — I > 1. 



(48) 



(49) 



Equation (49) implies that the position of the microcan- 
tilever can be very sensitive to the incoming ac signals 
if their frequency is close to one of the plasma resonant 
frequencies and the quality factor of the plasma oscilla- 
tions Q p oc il p /v ^> 1. Hence the micromashined HEMT 
under consideration can serve as a mechanical resonant 
detector of microwave and terahertz radiation. 

One may expect that at sufficiently strong ac signals, 
the microcantilever can be pulled-in to the surface of the 
isolating layer. Assuming, for simplicity that Vb = 0, 
we obtain the following condition of the microcantilever 
pull-in under the effect of the ac voltage: 



6 ZY> 



VoJ - 27^ ' 
Using the estimate for maxJ^, we obtain 



(50) 



so that the minimum ac pull-in voltage can be estimated 
as 



min 5v {pull ~ m) ~ ( -J- V - 



a, 



(52) 



Equation (54) shows that when the quality factor of the 
plasma oscillations is large, the microcantilever pull-in 
might occur at fairly modest ac signals. However, one 
may assume that the real situation is more complex be- 
cause the transition of the microcantilever to the position 
corresponding to its pull-in to the isolating layer should 
be accompanied by the channel depletion and significant 
change in the resonant plasma frequency. Due to this, the 
dynamic of the microcantilever pull-in out of the scope 
of this paper. 



VII. CONCLUSIONS 

We developed a model for a micromachincd HEMT 
with the microcantilever serving as the gate. The model 
is based on an equation of mechanical motion of the mi- 
crocantilever in the point-mass approximation accompa- 
nied by hydrodynamic equations describing distributed 
electron plasma phenomena in the 2DEG channel. Us- 
ing this model, we calculated the response function de- 
scribing the amplitude microcantilever oscillations and 
the output electric signal as functions of the signal fre- 
quency and the bias voltage for the devices with different 
parameters. We found the voltage dependences of the 
frequency of the mechanical resonance and its damping. 
It was demonstrated that the amplitudes of the mechan- 
ical oscillations and output electric signal exhibit pro- 
nounced maxima at the voltages close to the voltage of 
the 2DEG channel depletion. However, further increase 
in the bias voltage results in a drastic drop of the me- 
chanical and electrical response. We showed also that 
at the frequency corresponding to the plasma resonance 
the ac electric field between the microcantilever and the 
2DEG channel can be rather strong. This can result in 
a significant variation of the microcantilever position by 
incoming high-frequency (terahertz) signals. 

This work was partially supported by the Grant-in- 
Aid for Scientific Research (S) from the Japan Society 
for Promotion of Science, Japan. The work at RPI was 
partially supported by the Office of Naval Research, USA. 



Appendix A. Stationary states and pull-in and 
depletion voltages 

When the voltage applied between the gate and side 
contacts is constant, eqs. (2) - (5) yield u = u = 0, 
ip = fo = 0, and 



So 



v -v s 



4ire[Z a - w(l - ac" 1 )] 
with Zq governed by the following equation: 



Z a = W 



L g D 



V 2 



AirnlMlZo-wil-as- 1 ) 



(Al) 



(A2) 
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Introducing 

%= ^2nnlMW(W 2 /L g D), 
eq. (2) can be presented in the following form: 



W 



1 



2 \V& 



W 



Zo-wil-x- 1 ) 



(A3) 



The states with Zq satisfying eqs. (A2) and (A3) exist if 
\Vo\ <V { 




«,/ _1 
W V ac 



In 




(A4) 



is the so-called pull-in voltage [l| . One can find that when 

\Vq\ > Vq PuLI m , eq. (A3), does not have roots. In this 
case, the only existing (and stable) stationary state cor- 
responds to the attachment of the microcantilcver to the 

(pull— in) 



At |V | = V ( 



o 



isolation solid layer, i.e., to Zq = w. 
the microcantilever position is Zq = Zq PuLI in ' , where 



Z, 



(pull — in) 


W 



1 w 
3W 



(A5) 



When I Vol < V 



(pull — in) 







one with Zq < Z { 



?(pull—in) 
J 



(pull— in) 




eq. (A3) has two solutions: 
(unstable) and one with 



< Z < W (stable). At |V | < V [ 



(pull— in) 



the position of the microcantilever in the stable state is 
given by 



Zo _ 1 _l(Vb 
W~ 2\% 



(A6) 



There exists the stable state (not governed by eq. (A3)) 
with Zq = w as well . 

As follows from eqs. (A2) - (A3), Zq decreases with in- 
creasing | Vb| disregarding the sign of the gate voltage Vq. 
However, at Vq < 0, the 2DEG channel can be fully de- 
pleted, sothatEo = 0. The depletion voltage (or HEMT's 

threshold voltage) Vq P1 ' corresponding to £ = 0, as 
can be found from eqs. (Al) and (A3), is given approxi- 
mately by 



(depl) 



This voltage corresponds to 



-4ireE d W. 



7 (depl) 

'a 

W 



Vq 



V 



(depl) 



(A7) 



(A8) 



L g D 



25 /jhi 2 , W = 0.5 (im, M = 
= 1-10 MHz, we obtain V^ 



V 



(depl) 



9 - 90 V, and Vj 



(pull— in) 



18 



10" io g, and 
■ 33 - 330 V, 
- 180 V. Since 



Vq and, hence, Vq 



(pull— in) 



strongly decrease with decreas- 



ing W, their values can be markedly smaller than those 
obtained in the above estimates. The ratio oi\V r 



(depl) I 



to 



V 



(pull— in) 







can be presented as 



\V 



(depl) I 



(pull — in) 



< 2lTxe 2 Y? d L g D 

nlMw 



(A9) 



When Vq < and |Vb| > Vq , the charge densities 
in the microcantilever and channel do not change with 
varying Vq] they are equal to S^. In this case, eq. (A3) 
should be replaced by 



Zo 
W 



1 (Vr 



(depl) 



V 



Vo 



w 



Zo - W(l -83" 1 ) 



(A10) 



As a result, for the pull-in voltage under the condition of 
the channel depletion we obtain 



V, 



pull — in 



1 Vq 



2y r 



(depl) 



W 

w 



1 Vo 



*Vr 



and (compare with eq. (A5)) 



7 (pull—in) 

•a 

W 



w 

w 



1 - 



(depl) 

(All) 
(A12) 



Appendix B. Voltage and mechanical control of the 
plasma resonances 

The expression for the characteristic plasma frequency 
il p of the gated 2DEG channel given by eq. (20) is some- 
what different from that obtained previously for the stan- 
dard HEMTs. This frequency in micromachincd HEMTs 
exhibits different voltage dependence. As follows from 
eq. (20), the resonant plasma frequency f2 p depends not 
only on the electron density in the 2DEG channel but also 
on the spacing, Zq, between the microcantilever (gate) 
and 2DEG. This opens up the possibility of a mechani- 
cal control of the plasma frequency and, hence, different 
characteristics of the pertinent terahertz devices. 

One can obtain the following dependence of the char- 
acteristic plasma frequency on the microcantilever dis- 
placement and the bias voltage. 




(Bl) 



The second terms in the right-hand side of Eq. (A8) is 
rather small. Indeed, if = (1 — 10) x 10 11 cm~ 2 , 



If Vq is fixed, the dependence of the characteristic plasma 
frequency on the microcantilever displacement is given by 
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eq. (41): fl p oc \[Z§. Taking into account that the micro- 
cantilcvcr displacement depends on the bias voltage (see 
eq. (A6) from Appendix A), we find the following voltage 
dependence of the characteristic plasma frequency: 

^^V i+ KtH(t) 2 ' (B2) 



where (3 = Vo/\Vq ep \. It differs from the pertinent 
formula for the standard HEMTs by the third term under 
the square root in the right-hand side of eq. (B2). 
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